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^ . ABSTRACT 

A\ms. We present the first multi-epoch A'-band spectro-interferometric observations of the carbon-rich Mira variable V Oph using MIDI at the 
ESO's Very Large Telescope Interferometer. Our aim is to study temporal variations of physical properties of the outer atmosphere and the 
circumstellar dust shell based on spectrally-dispersed W-band visibilities over the C2H2 (+HCN) features and the dust emission. 
, Methods. Our MIDI observations were carried out at three different phases 0.18, 0.49, and 0.65, with three different baselines (projected 
C ■ baseline lengths of 42-124 m) using four 8.2 m Unit Telescopes (UT2-UT4, UTI-UT4, and UT2-UT3 baseline configurations), 
^p " Results. The wavelength dependence of the uniform-disk diameters obtained at all epochs is characterized by a roughly constant region between 
8 and 10 yum with a slight dip centered at ~9.5 yum and a gradual increase longward of 10 yum. These A'-band angular sizes are significantly 
O I larger than the estimated photospheric size of V Oph. The angular sizes observed at different epochs reveal that the object appears smaller at 
' ^ phase 0.49 (minimum light) with uniform-disk diameters of ~5-I2 mas than at phases 0.18 (~I2-20 mas) and 0.65 (~9-15 mas). We interpret 
^ . these results with a model consisting of optically thick C2H2 layers and an optically thin dust shell. Our modeling suggests that the C2H2 layers 
■ around V Oph are more extended (~ 1.7-1.8 R^,) at phases 0.18 and 0.65 than at phase 0.49 (~1.4 R^,) and that the C2H2 column densities appear 
to be the smallest at phase 0.49. We also find that the dust shell consists of amorphous carbon and SiC with an inner radius of ~2.5 and the 
• • total optical depths of Ty ~ 0.6-0.9 (Tn.3,jm ~ 0.003 and 0.004 for amorphous carbon and SiC, respectively) found at phases 0. 18 and 0.65 are 
. 1^ higher than the value obtained at phase 0.49, Tv~0.3 (th s^n, r; 0.001 and 0.002 for amorphous carbon and SiC, respectively). 

Conclusions. Our MIDI observations and modeling indicate that carbon-rich Miras also have extended layers of polyatomic molecules as 
» I previously confirmed in oxygen-rich Miras. The temporal variation of the A'-band angular size is largely governed by the variations of the 
^ opacity and the geometrical extension of the C2H2 layers and the dust shell, and consequently, this masks the size variation of the photosphere. 
Also, the observed weakness of the mid-infrared C2II2 absorption in carbon-rich Miras can be explained by the emission from the extended 
C2H2 layers and the dust shell. 

Key words, infrared: stars - techniques: interferometric - stars: circumstellar matter - stars: carbon - stars: AGB and post-AGB - stars: 
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1. Introduction 

The driving mechanism of mass outflows in Mira variables 
has not yet been fuUy understood. Recent progress in optical 
and infrared interferometric techniques has been contributing 
to studies of the region between the top of the photosphere 
and the innermost region of the circumstellar dust sheU, exactly 
where mass outflows are expected to be initiated. Near-infrared 
interferometric observations of several oxygen-rich and S-type 

Send offprint requests to: K. Ohnaka 

* Based on observations made with the Very Large Telescope 
Interferometer of the European Southern Observatory. Program ID: 
075.D-0607 

** Table. 4 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A 



(C/O » 1) Mira variables (Mennesson et al. 120021 Perrin et 
al. 120041 Ireland et al. |20D4l Woodruff' et al. |2004l Fedele 
et al. '2005) have turned out to be consistent with the pres- 
ence of dense molecular layers extending to ~2 i?*, which 
were introduced to explain spectroscopic observations (e.g., 
Hinkle et al. [T9791 Tsuji et al. [19971 Yamamura et al. \T^99\ . 
Furthermore, based on semi-empirical models, Ohnaka ( |2004al l 
quantitatively shows that near-infrared (K and L' bands) and 
mid-infrared (11 yum) interferometric and spectroscopic obser- 
vations of three oxygen-rich and S-type Miras (o Cet, R Leo, 
and;^' Cyg) can be explained by dense warm water vapor layers 
extending to -2.5 with temperatures of 1200-2000 K and 
H2O column densities of the order of 10~' cm"^. These dense 
water vapor layers are responsible for the increase of the an- 
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gular size from the near-infrared to the mid-infrared observed 
toward oxygen-rich Miras. 

However, these studies have been Hmited to oxygen-rich or 
S-type Mira variables up to now, and the physical properties of 
the outer atmosphere of carbon stars and their temporal varia- 
tions have not yet been well probed. Observations of carbon- 
rich Miras as well as non-Mira carbon stars using the Short 
Wavelength Spectrometer (SWS) onboard the Infrared Space 
Observatory (ISO) revealed that while the spectra of carbon 
stars shortward of ~5 /im can be fairly reproduced by non-gray 
hydrostatic or dynamical model atmospheres, these models en- 
tirely fail to explain the observed spectra longward of ~5 fj.m, as 
demonstrated by J0rgensen et al. (|20001 l and Gautchy-Loidl et 
al. ( I2004I I. These models predict very strong absorption due to 
C2H2 and HCN at 7 and 14 fim, but the observed ISO spectra of 
carbon stars show only weak absorption due to these molecular 
species. As a possible solution for this discrepancy, Gautchy- 
Loidl et al. (2004) show that dynamical models with dust- 
driven stellar winds can provide a qualitative explanation of the 
absence of strong absorption features due to C2H2 and HCN at 
least for stars with mass loss rates of ~10"'' Mq yr"', although 
these models cannot yet reproduce the observed strengths of 
the molecular features quantitatively. On the other hand, Aoki 
et al. (119981 11999I I propose that emission from extended warm 
molecular layers containing C2H2 and HCN may be responsi- 
ble for the discrepancy between the ISO spectra and the photo- 
spheric models. 

Infrared interferometry provides us with an excellent op- 
portunity to probe the nature of this controversial region — 
the outer atmosphere and the innermost region of the dust 
shell of carbon stars. Van Belle et al. ( 119971 ) observed two 
carbon-rich Miras U Cyg and V Cyg at 2-3 epochs in the K 
band using the Infrared Optical Telescope Array (IOTA), but 
the accuracy of the derived angular sizes as well as the pres- 
ence of molecular absorption, which affects the apparent stel- 
lar size in these broad-band observations, makes the deduc- 
tion of temporal variations difficult. Thompson et al. ( 120021 ) 
brought about a clearer picture of the phase dependence of 
the angular size of the carbon-rich Mira RZ Peg based on 
multi-epoch narrow-band observations in the TiT-band using the 
Palomar Testbed Interferometer (PTI): the uniform-disk diam- 
eters measured over two pulsational cycles show a variation 
from 2.3 to 3.0 mas. Their measurements also show that the 
observed uniform-disk diameter monotonically increases from 
2.0 to 2.4 fj.m, which can be interpreted as that the strong CO 
first overtone bands longward of 2.3 fim make the star appear 
larger than at shorter wavelengths where only relatively weak 
(but crowded) CN lines are present. 

However, mid-infrared interferometric observations of car- 
bon stars (particularly, not very dusty ones) are still very 
scarce, let alone multi-epoch measurements to follow tem- 
poral variations of the physical properties of the photo- 
sphere and the circumstellar environment. The MID-infrared 
Interferometric instrument (MIDI) at the ESO's Very Large 
Telescope Interferometer (VLTI) is well suited for a study of 
the circumstellar environment close to the star It enables us to 
obtain visibilities with a spectral resolution of 30 or 230 from 
8 to 13 jum, where the absorption features due to C2H2 and 



HCN as well as the SiC and (featureless) amorphous carbon 
dust emission are observed for carbon stars. 

In this paper, we present the first multi-epoch A^-band 
spectro-interferometric observations of the carbon-rich Mira 
star V Oph using MIDI. Our object V Oph has a period of 
297 days (Samus et al. I2004 i and the V magnitude varies 
from 7.3 to 11.6 (Fn/im ~ 20 Jy). Its spectral energy distri- 
bution (SED) suggests that the circumstellar dust shell is op- 
tically thin with quite a low mass loss rate of ~10"^ Mq yr"' 
(Groenewegen et al. 1998). The absence of a thick dust shell 
makes V Oph suitable for probing the presence or absence of 
extended molecular layers in carbon-rich Miras. 

2. MIDI observations 

V Oph was observed with MIDI on six nights between April 
and September 2005 in service mode (Program ID: 075-D0607, 
P.I.: K. Ohnaka), as summarized in Table [T] A prism with a 
spectral resolution of A/AA ^ 30 at 10 jum was used to obtain 
spectrally dispersed fringes between 8 and 13 yum. A detailed 
description of the observing procedure is given in Przygodda 
et al. (2003j), Leinert et al. (120041 ). and Chesneau et al. (I2005al 
I2005bl) . 

We use two different MIDI data reduction packages: MIA 
developed at the Max-Planck-Institut fiir Astronomic and EWS 
developed at the Leiden Observatory (MIA-hEWS, ver.l.4il|. 
While the MIA package is based on the power spectrum 
analysis, which measures the total power of observed fringes 
(Leinert et al. 120041) . the EWS software first corrects for optical 
path differences (instrumental as well as atmospheric delays) 
in each scan, and then, the fringes are coherently added (Jaffe 
120041) . We derive the interferometer transfer function at each 
wavelength (i.e., at each spectral channel) between 8 and 13 fim 
by observing calibrators whose angular diameters are known 
(see Table |2|. In order to calibrate the visibility of V Oph at 
a given wavelength, we use the mean of the transfer function 
values derived from all the calibrators observed on the same 
night as V Oph, and the errors of the calibrated visibilities are 
estimated from the standard deviation of the transfer function 
values at each wavelength. The results obtained with the MIA 
and EWS packages show good agreement, and we only present 
the result reduced with the EWS package below. The relative 
errors of the calibrated visibilities are typically +10-15% (Icr), 
and the error sources are described in Ohnaka et al. (120051) . 

Figure[T]shows the calibrated visibilities of V Oph obtained 
on six nights (the complete set of the observed visibilities are 
electronically available in Table 4). The observed A^-band visi- 
bilities except for the data #4 taken with the shortest baseline on 
2005 July 23 are characterized by a monotonic increase from 
8 to ~9.5 fj.m and a gradual decrease longward of ~10 yum. 
Figures |2h-c show the observed visibilities plotted as a func- 
tion of spatial frequency at three representative wavelengths be- 
tween 8 and 13 yum. In the modeling described below, the data 
sets are binned into three epochs: epoch 1 (2005 Apr. 17 + 18, 
phase 0.18), epoch 2 (2005 Jul. 19 + 23, phase 0.49), and epoch 

' Available at http://www.mpia-hd.mpg.de/MID ISOFT/| and 
|http ://vy vyw. strw.leidenuniv. nl/' nevec/MIDI/index. html 
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B =B5.6in, P.A. = 85.4° (UT2-UT4) 



2005 Apr. 18, phase = 0,18 
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2005 Jul. 19, phase = 0.49 
fl„=124.2m, P. A. = 65. 6° {UT1-UT4) 



2005 Jul. 23, phase = 0.49 

B =41.9m, P.A. = 25.7° {UT2-UT3) 
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2005 Aug. 24, phase = 0.61 

=87.6m, P.A. = 84.3° (UT2-UT4) 



9 10 11 12 13 14 

A (/Am) 



2005 Sep. 18, phase = 0.69 

B =120.0m, P.A.=66.0° (UT1-UT4) 



9 10 11 12 13 14 

A (/xm) 



Fig. 1. A'-band visibilities of V Oph ob- 
served with MIDI on six nights between 
April and September 2005. 



Table 1. Summary of MIDI observations of V Oph: night, time 
of observation (UTC=Coordinated Universal Time), airmass 
(AM), telescope configurations (Tel., UT=Unit Telescope), 
projected baseline length Bp, position angle of the baseline vec- 
tor on the sky (P.A.), and variability phase estimated from the 
visual light curve compiled by the American Association of 
Variable Star Observers (AAVSO). 



# 


Night 


^obs 


AM 


Tel. 


5p 


P.A. 


phase 




(2005) 


(UTC) 






(m) 


n 




1 


Apr. 17 


08:39:02 


1.1 


UT2-4 


85.6 


85.4 


0.18 


2 


Apr. 18 


04:51:38 


1.3 


UT2-4 


71.8 


70.1 


0.18 


3 


Jul. 19 


03:03:08 


1.1 


UTl-4 


124.2 


65.6 


0.49 


4 


Jul. 23 


22:58:21 


1.2 


UT2-3 


41.9 


25.7 


0.49 


5 


Aug. 24 


23:49:30 


1.1 


UT2-4 


87.6 


84.3 


0.61 


6 


Sep. 18 


23:27:32 


1.2 


UTI-4 


120.0 


66.0 


0.69 



3 (2005 Aug. 24 + Sep. 18, phase 0.61-0.69 with an average 
of 0.65). The data sets #1 and #5, which were obtained four 
months apart but at almost the same baseline lengths and posi- 
tion angles, reveal a temporal variation of visibility, particularly 
at 8.3 jum. However, the errors of the measured visibilities be- 



come larger at longer wavelengths, which makes the detection 
of temporal variations less conclusive. For example, the large 
errors of the visibilities measured at 12.5 fim hamper us from 
definitively detecting a temporal variation. On the other hand, 
the data sets #3 and #6, which were obtained two months apart 
at approximately the same baseline lengths and position angles, 
clearly indicate a temporal variation not only at 8.3 fim but also 
at 10 and 12.5 yum. 

At each epoch, we fit the visibilities at a given wavelength 
with a uniform disk as shown in Figs. |2^-c. While the differ- 
ences in position angle for the data sets of the first and third 
epochs are not very large (they differ by 15.3° and 18.3°, re- 
spectively), the data sets of the second epoch (#3 and #4) were 
obtained with quite different position angles of 65.6° and 25.7°. 
We tentatively fit the visibility points of these data sets with a 
single uniform disk, regardless the position angles, because the 
difference in baseline length between two data sets makes it 
impossible to draw a conclusion about the possible presence 
or absence of asymmetries. The resulting uniform-disk diame- 
ters derived for three epochs are plotted as a function of wave- 
length in Fig. |2]l, which illustrates a temporal variation of the 
angular size of V Oph: the object appears smaller at minimum 
light (phase 0.49) than at post-maximum (phase 0.18) or post- 
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Fig. 2. A^-band visibilities of V Oph plotted as a function of spa- 
tial frequency at three representative wavelengths, a: 8.3 /im, 
b: 10.0 fim, and c: 12.5 fim. In these panels, the filled circles, 
diamonds, and triangles represent the data obtained at epoch 1 
(phase 0.18), epoch 2 (phase 0.49), and epoch 3 (phase 0.65), 
respectively, as defined in Sect. |2] The fit with uniform disks is 
also plotted for three epochs. The data set number of each vis- 
ibility point as given in Table [1] is indicated, d: Uniform-disk 
(UD) diameters derived by fitting the visibility points obtained 
at each epoch are plotted as a function of wavelength. The filled 
circles, diamonds, and triangles represent the data obtained at 
epoch 1, 2, and 3, respectively. 



Table 2. List of calibrators used for the calibration of the 
V Oph data, together with 12 fim fluxes (F12), uniform-disk 
diameters (liuo)- the date as well as the time stamp (fobs) of 
MIDI observations, and airmass (AM). The uniform-disk di- 
ameters were taken from the CalVin list available at ESO 
(http://www.eso.org/observing/etc/). The data sets used for 
spectrophotometric caUbration of the spectra of V Oph are 
marked with f. 



Calibrator 






Night 


'obs 


AM 




(Jy) 


(mas) 


(2005) 


(UTC) 




HD 120404 


13.3 


3.03 ± 0.24 


Apr. 17 


02:48:03 


1.5 


HD 142804 


11.7 


2.80 ± 0.08 


Apr. 17 


04:21:08 


1.2 








Apr. 18 


03:26:25 


1.5 


HD 169767 


8.8 


2.16±0.11 


Apr. 17 


06:08:59 


1.4 








Apr. 17 


07:18:00 


1.2 








Apr. 17 


07:53:19 


1.2 


HD150052 


8.9 


2.40 + 0.07 


Apr. 17 


08:17:33+ 


1.1 








Apr. 18 


04:29:57 


1.4 








Jul. 19 


02:43:38"'" 


1.1 








Jul.23 


23:20:17 


1.2 








Aug.24 


23:22:48+ 


1.0 








Sep. 18 


23:51:54+ 


1.3 


HD188512 


9.0 


2.07 + 0.10 


Apr. 17 


09:04:41 


1.3 








Jul.23 


02:58:55 


1.3 








Jul.23 


04:03:14 


1.2 








Jul.23 


06:37:58 


1.4 


HD169916 


31.2 


3.90 + 0.21 


Apr. 17 


10:08:16"+ 


1.0 








Apr. 18 


07:24:08 


1.1 








Jul. 19 


05:46:15 


1.2 








Jul.23 


00:04:56+ 


1.3 








Jul.23 


05:49:35 


1.3 


HD37160 


9.4 


2.08 ± 0.20 


Apr. 18 


23:26:37 


1.7 


HD 122451 


11.4 


0.96 + 0.07 


Apr. 18 


02:26:48 


1.4 


HD152885 


12.6 


2.88 ± 0.09 


Apr. 18 


05:12:20 


1.2 








Jul.23 


00:48:11 


1.0 


HD 157236 


10.5 


2.58 ±0.14 


Apr. 18 


06:21:08+ 


1.1 








Apr. 18 


08:16:35 


1.0 








Aug.24 


01:40:04 


1.1 


HD 160668 


8.4 


2.30 ±0.10 


Apr. 18 


06:43:00 


1.1 


HD168723 


17.3 


2.88 ±0.13 


Apr. 18 


09:17:26 


1.1 


HD152820 


10.9 


2.63 ±0.16 


Jul. 19 


04:14:50 


1.2 








Jul.23 


02:15:18 


1.0 








Jul.23 


05:01:41 


1.5 








Sep. 18 


00:40:11 + 


1.3 


HD178345 


12.5 


2.50 ± 0.07 


Aug.24 


04:24:35 


1.3 


HD165135 


23.4 


3.38 ±0.16 


Sep. 18 


01:30:20 


1.2 



minimum (phase 0.65). It should be stressed, however, that we 
use uniform-disk fits to obtain some kind of representative an- 
gular size of the object and that it can be very misleading to 
deduce quantitatively temporal variations from these uniform- 
disk diameters, particularly if the real intensity distribution of 
the object is very different from a uniform disk. 

The photospheric angular size of V Oph (i.e., the stellar 
angular diameter measured in the continuum) has not been di- 
rectly measured, but we can estimate it from the bolometric 
flux and effective temperature, if the radiation is assumed to be 
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Fig. 3. A^-band spectra of V Oph obtained at three epochs. The 
filled circles, diamonds, and triangles represent the spectra ob- 
tained at epochs 1, 2, and 3, respectively. The shaded region 
represents the wavelength range severely aff'ected by the tel- 
luric ozone absorption centered at 9.6 /im. 



isotropic. The bolometric flux of V Oph near maximum light is 
derived to be 4.4 x 10"'" W m"^ from photometric data avail- 
able in the literature, while the effective temperature is esti- 
mated to be 2600-3000 K (see Sect. [3] for discussion on pho- 
tometric data and effective temperature). These bolometric flux 
and efl'ective temperatures translate into a photospheric angu- 
lar size of 4-6 mas. It means that the A^-band angular sizes 
obtained at three epochs are significantly larger than the pho- 
tospheric size in most cases. Only the 9-10 //m uniform-disk 
diameters obtained at phase 0.49 are rather close to the photo- 
spheric size. Therefore, for the carbon-rich Mira V Oph, we de- 
tect the same trend as found in oxygen-rich Mira stars: the mid- 
infrared angular sizes are much larger than the photospheric 
angular size. 

We also extract the absolutely calibrated A^-band spectra of 
V Oph from the MIDI data by dividing the observed spectra 
of V Oph by that of a spectrophotometric standard star and 
then multiplying the absolutely calibrated spectrum of this lat- 
ter standard star. Interferometric calibrators observed close in 
time to V Oph at similar airmasses are used for extracting the 
spectra of V Oph as marked with daggers in Table|2l The abso- 
lutely calibrated spectra of these spectrophotometric standard 
stars are taken from Cohen et al. ( 1999i ) if available. Otherwise, 
the absolutely calibrated spectrum of a standard star is calcu- 
lated as follows. First, we select a star in Cohen et al. (119991) 
with the same spectral type and luminosity class as our stan- 
dard star Then the flux level of the spectrum of this "template" 
star is scaled using its angular diameter given in Cohen et al. 
( 11999 1) and that of our standard star given in Table |2l and the 
resulting spectrum is used for spectrophotometric calibration. 
Finally, we average all the resulting spectra obtained at each 
epoch. The absolutely calibrated spectra of V Oph observed at 
three epochs are shown in Fig. [3] The errors of the MIDI spec- 
tra are difficult to estimate properly from the present data sets, 
because the number of standard stars observed close in time 
to V Oph and at similar airmasses is quite limited (3, 2, and 
3 standard stars at epoch 1, 2, and 3, respectively). However, 
we estimate the errors of the absolutely calibrated spectra to be 
10-20%, based on other studies which derived A^-band spectra 



from MIDI data (e.g., Kervella et al. |2006t Poncelet et al. l2006l 
Quanz etal.|200g). 

Figure [3] reveals that the spectra are characterized by ab- 
sorption between 8 and 9 jim and a broad emission feature cen- 
tered at ~ 11 .3 jum. The feature seen between 8 and 9 fim is actu- 
ally the tail of the broad absorption feature centered at ~7.5 fim 
due to the C2H2 V4-1-V5 band and the HCN 2v2 band. However, 
as shown in Sect. 14.11 the HCN 2v2 band extends only up to 
~8.3 yum, which leaves C2H2 as the dominant opacity source 
between 8 and 9 fim. The broad emission feature at ~11.3 jum 
is attributed to SiC dust. In addition to these conspicuous fea- 
tures within the A^ band, a broad absorption feature centered 
at ~14 fim due to the C2H2 V5 and HCN V2 bands is identi- 
fied in ISO spectra of other carbon-rich Mira variables (e.g., 
Yamamura et al. [T997l Hron et al. 1998; Aoki et al.fT999^. This 
absorption feature extends down to ~1 1 fim, where it becomes 
blended with the SiC emission feature. Comparison between 
the spectra at phases 0.18 and 0.49 reveals that the 8 //m ab- 
sorption feature becomes somewhat stronger at minimum light 
than at post-maximum. This trend is similar to that found in 
the Mira-like carbon star R Scl by Hron et al. (|19981 l. We will 
discuss this variation of the 8 fim feature in Sect. 14.21 

3. Modeling with a dust shell 

As briefly discussed in Sect. |2] the A^-band uniform-disk diam- 
eters measured with MIDI are remarkably larger than the esti- 
mated photospheric angular size, which suggests the contribu- 
tion of circumstellar material in the mid-infrared angular size. 
This contribution can be of the molecular and/or dust origin. 
We first attempt to interpret the observed A^-band visibilities 
and spectra with dust shell models. 

We calculate spherically symmetric dust shell models, 
using the multi-dimensional Monte Carlo code mcsim_mpi 
(Ohnaka et al. 2006a), which computes dust temperatures in 
radiative equilibrium and generates model SEDs. Then the 
model image at each wavelength within the A^ band is pro- 
duced with the ray-trace method as described in Ohnaka et 
al. (|2006al l, and the corresponding monochromatic visibility is 
calculated by taking the Hankel transform of the model image. 
The monochromatic visibilities are then spectrally convolved to 
a resolution of 30 to match the MIDI observations. The effec- 
tive temperature of V Oph derived by various authors ranges 
from 2650 to 3010 K (Ohnaka & Tsuji [T9961 Bergeat et al. 
120011 Groenewegen2Q06), but the phase dependence of effec- 
tive temperature of V Oph is not well established. Therefore, in 
the dust shell modeling below, we adopt a fixed effective tem- 
perature of 2800 K, which is roughly the mean of the above 
range, and the central star is approximated with the blackbody 
of 2800 K. We assume that the dust shell consists of amorphous 
carbon and SiC, following the results of Groenewegen et al. 
( 119981 ) and Groenewegen (120061) . who concluded that the dust 
chemical composition around V Oph is dominated by amor- 
phous carbon (85%) with a minor contribution of SiC (15%). 
We compute a grid of models with the following parameter 
ranges: n„ (R^) = 2.5 ... 5.0 (Ar;,, = 0.5), t™'= = 0.2 ... 0.8 
(At™'^ = 0.2), where ri„ and r™"^ denote the inner boundary 
radius (assumed to be the same for amorphous carbon and SiC) 
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Fig. 4. Dust shell model containing amor- 
phous carbon and SiC which well repro- 
duces the observed SED for phases near 
maximum light. The dust shell parameters 
are = 2.5 R^,, = 0.6, and t^"= = 0.26. 
a: Comparison of the observed and model 
SEDs. The observed SEDs are represented 
with the filled diamonds, while the model 
SED is represented with the solid line, b: 
Comparison of the observed and model A'- 
band spectra. The filled circles represent 
the MIDI spectrum obtained at phase 0.18, 
while the solid line represents the model 
prediction, c: The filled diamonds and trian- 
gles represent the A'-band visibilities mea- 
sured with projected baseline lengths of 
85.6 and 71.8 m, respectively, while the 
solid and dashed lines represent the corre- 
sponding model visibilities, d-f: The solid 
lines represent the model visibilities at 8.3, 
10.0, and 12.5 /jm plotted as a function 
of spatial frequency. The observed data are 
plotted with the filled circles. 



and the optical depth of amorphous carbon dust at 0.55 /jm, re- 
spectively. The relative contribution of SiC, t^''--/(t™'^ + t^'~) 
is varied from 0.1 to 0.4 with an increment of 0.1, where t^'~ 
denotes the optical depth of SiC dust at 0.55 fim. The opaci- 
ties of amorphous carbon and SiC are calculated from complex 
refractive indices presented by Rouleau & Martin (|1991i ACl 
sample) and by Pegourie (1988 ), respectively, assuming a sin- 
gle grain size of 0.1 /jm for both species. The grain density 
distribution is assumed to be proportional to r^^, and the outer 
radius of the shell is fixed to 10^ x ^n. 

Figure H] shows a comparison between the observational 
data and a model which can well reproduce the SED ob- 
served near maximum light. We derived the observed SED 
from optical and infrared photometric data available in the 
literature (Landolt fT968l 119731, Walker IT9791 Noguchi et al. 
[19811 Whitelock et al. 2000; IRAS Point Source Catalog) and 
dereddened with A v = 0.6 (Groenewegen et al. |1998) using the 
method of Savage & Mathis (fT979] l with Ay = 3.1E{B - V). 
Unfortunately, most of the photometric data, particularly in 
the optical, were obtained at phases near maximum light, and 
therefore, we can compare model SEDs with the observed one 
only for phase 0.18. As can be seen in Figs. andl^j, the 
SED at maximum light as well as the MIDI spectrum ob- 
tained at phase 0.18 can be explained by the model with = 
2.5 R^, 4™ = 0.6 (0.0038 at 1 1.3 fim), and 4"= = 0.26 (0.0078 
at 11.3 fim). The temperatures of amorphous carbon and SiC 
dust at the inner boundary are ~1580 K, which roughly corre- 
sponds to the condensation temperature of amorphous carbon 
predicted by dynamical models (e.g., Nowotny et al. I2005al l 
and is in agreement with the value that Groenewegen et al. 
d 19981 1 derived from a SED-fitting (note, however, that conden- 
sation temperatures generally depend on detailed grain nucle- 
ation and growth processes). The model cannot reproduce the 
8 jum feature, because the star is approximated with a black- 



body without a detailed model atmosphere. This is also the rea- 
son the model overestimates the flux in the optical where strong 
molecular absorption is present. Comparison of the visibilities 
shown in Fig.|4j; reveals that the dust shell model can roughly 
reproduce the observed visibility levels. The dip at ~11.3 fim 
seen in the model visibilities corresponds to the position of the 
SiC emission feature. The flux contribution of the dust shell 
becomes larger due to the SiC emission feature at these wave- 
lengths, which leads to lower visibilities. It should be noted that 
the scaling factor between a model SED and the observed one is 
given by {Ri^jdf-, where d is the distance of V Oph. The model 
visibilities are scaled with the angular radius of the central star 
(Ri,/d) derived by fitting the model SEDs to the observed data. 

However, there are noticeable discrepancies between the 
model visibilities and the MIDI data. Firstly, the models can- 
not reproduce the observed wavelength dependence of visibil- 
ity between 8 and 9 /im as well as long ward of 12 yum, where 
the C2H2 (and HCN) absorption features are present, although 
the error bars of the observed visibilities are rather large in the 
latter wavelength range. The fact that the observed visibilities 
are smaller than the model at 8-9 yum and ^12 jjm suggests 
that the observed angular sizes are larger than that predicted by 
the dust shell model. Secondly, the models cannot reproduce 
the shape of the observed visibilities as a function of spatial 
frequency, which can be seen in Fig. |4jl: the model visibility 
at 8.3 jjm is somewhat too flat compared to those observed. 
Since the dust shell is over-resolved with the baselines used in 
our MIDI observations (the dust shell visibility component is 
seen as steep visibility drops at spatial frequencies from to 
-5 cycles/arcsec in Figs.|4}l-f), the visibility shapes measured 
by our observations correspond to the central star component. 
Therefore, this disagreement of visibility shape means that the 
observed A^-band angular size of the central star is larger than 
that used in the model. However, the angular size of the cen- 
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tral star is derived by fitting the model SED to the observed 
one, and therefore, a mere significant change of the angular 
size of the central star would lead to a poorer match to the ob- 
served SED including the A^-band spectrum. Changes in effec- 
tive temperature, dust chemistry, and density gradient cannot 
explain the observed visibility shapes, either. Therefore, it can 
be postulated that there is some additional component close to 
the star, which makes the central source appear larger than the 
star itself. 

4. Modeling with C2I-I2 layers and a dust shell 

As in oxygen-rich Mira variables which have dense H2O layers 
extending to 2-3 R*, it is possible that similar layers consisting 
of polyatomic molecules exist around carbon-rich Miras, which 
makes the central star appear larger than its photospheric size. 
In carbon-rich atmospheres, C2H2 and HCN are the most abun- 
dant polyatomic species with spectral features in the band, 
which makes them the most likely candidates responsible for 
the larger apparent size of the central star component. In this 
section, we attempt to explain the observed data with models 
consisting of such molecular layers and an optically thin dust 
shell. 

4.1. Molecular data 

First, we describe molecular line opacity data necessary for 
computing models with C2H2 (+ HCN) layers. In the wave- 
length region between 8 and 9 ;um, the C2H2 V4-1-V5 band is 
the dominant opacity source. It is also necessary to take the 
C2H2 V5 band centered at 14 /im into account for modeling the 
A^-band visibilities and spectra, because this band extends to 
~1 1 yum. The HCN 2v2 band is centered at ~7 jim, and its tail 
extends to ~8.3 fim. The HCN V2 band, which is observed to- 
gether with the C2H2 V5 band in other carbon-rich Miras, is 
centered at 14 fim and extends to ~11 fim. Therefore, the con- 
tribution of HCN appears mainly at wavelengths between 1 1 
and 13 //m, blended with the C2H2 V5 band and the SiC dust 
feature. 

Despite its importance in various astrophysical environ- 
ments, the line list of C2H2 is far from complete, particularly 
for high temperatures relevant to stellar photospheres and outer 
atmospheres (temperatures of 1000-1600 K are derived in our 
modeling below). There are a number of laboratory measure- 
ments of the molecular constants of the V4H-V5 and V5 bands, but 
these experiments were performed at room temperatures. For 
example, Kabbadj et al. ( |1991| l measured the molecular con- 
stants up to V4 + Vs < 4, which means that only 10 and 6 
transitions among different (V4, V5) are included for the V5 and 
V4-1-V5 bands, respectively. Therefore, many hot bands expected 
in stellar and circumstellar environments are not included, and 
these lists are not appropriate for our modeling. On the other 
hand, Tsuji (1984) presents a totally different approach to deal 
with C2H2 opacity. He makes use of the band model method, 
which analytically calculates the absorption cross section at a 
given wavelength from basic molecular constants, instead of a 
detailed line list. Particularly, if the molecule is assumed to be 
a rigid-rotator-harmonic-oscillator (RRHO), calculations of ab- 
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Fig. 5. Emission spectra of C2H2 and HCN from a slab of 1000 
and 1500 K. The column densities of C2H2 and HCN are set 
to 10^° cm"^. The spectra are convolved with a spectral reso- 
lution of 30. The opacity of C2H2 is computed with the band 
model method as described in Sect. 14.11 while that of HCN is 
computed from the line list of Harris et al. (.2002 J . 



sorption cross sections can be easily performecd In the present 
work, we use this band model opacity computed in the RRHO 
approximation to represent the C2H2 opacity, because neither 
laboratory measurements nor ab initio calculations at high tem- 
peratures are yet available for C2H2 lines, and individual lines 
are not resolved with the MIDI'S spectral resolution of 30. 

Compared to C2H2, studies on HCN line lists including 
high vibrational levels are more advanced. Recently, Harris et 
al. ( 120021) have calculated a new, extensive ab initio line list of 
HCN/HNC including energy levels up to 18 000 cm We use 
their line list to calculate the HCN opacity in the modeling be- 
low. Figure |5] shows emission spectra of C2H2 and HCN from 
a slab of 1000 and 1500 K with the column densities of both 
molecules set to 10^° cm"^. The HCN opacity was calculated 
assuming local thermodynamical equilibrium (LTE) adopting 
a Gaussian profile with a FWHM of 5 km s ', and the opacity 
of C2H2 was calculated with the band model method with the 
RRHO approximation, which also implicitly assumes LTE. The 
spectra are convolved with a spectral resolution of 30 as used in 
the MIDI observations. The figure illustrates that C2H2 is the 
dominant opacity source at wavelengths between 8 and 9 jum 
as well as longward of 1 1 yum, while the contribution of HCN 
appears primarily at wavelengths longer than 1 1 /urn (and at the 
shorter edge of the band). Since this contribution of HCN is 
blended with the C2H2 V5 band and the SiC dust feature, the 
effects of HCN on spectra and visibilities tend to be masked 
by the C2H2 and SiC features, which hinders us from putting 
strong constraints on the HCN column density. Therefore, in 
the models presented below, we first consider C2H2 alone and 
examine the effects of HCN later 



^ Figure 3 of Tsuji ( I1984t . which shows the absorption cross sec- 
tion of the C2H2 V5 band, has an error in the numerical labels of the 
ordinate. The correct range of the ordinate is from -22 to -17 instead 
of -23 to -18. However, this is a mere graphical error, and the results 
of the calculations in Tsuji ([1984} were performed with the correct 
cross sections (Tsuji, priv. comm.). 
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Fig. 6. Comparison between observational data at phase 0.18 (post-maximum) and the best-fit model consisting of C2H2 layers 
and a dust shell. The model parameters are given in Table|3] a: The filled diamonds and triangles represent the A^-band visibilities 
measured at phase 0.18 with projected baseline lengths of 85.6 and 71.8 m, respectively, while the solid and dashed lines represent 
the corresponding model visibilities, b: The filled circles represent the uniform-disk diameters observed at phase 0.18, while the 
solid line represents the model prediction, which is derived by fitting the model visibilities predicted for baselines of 85.6 and 
71.8m with a uniform disk at each wavelength. The dotted line represents the photospheric angular diameter derived by fitting 
the model A^-band spectrum to the observed data, c: Comparison of the observed and model A^-band spectra. The filled circles 
represent the MIDI spectrum obtained at phase 0.18, while the solid line represents the model prediction, d-f: The solid lines 
represent the model visibilities at 8.3, 10.0, and 12.5 /urn plotted as a function of spatial frequency. The observed data are plotted 
with the filled circles, g: The C2H2 optical depths of the hot and cool layers (convolved to a spectral resolution of 30). h: The 
optical depths for amorphous carbon (AMC) and SiC dust (convolved to a spectral resolution of 30). 



4.2. Result of modeling with C2H2 layers and a dust 
shell 

We start from models consisting of one C2H2 layer with a con- 
stant temperature and density and an optically thin dust shell 
of amorphous carbon and SiC, similar to the model used by 
Ohnaka et al. ( |20051 l. In this model, the temperature, outer ra- 
dius, and the C2H2 column density of the layer are varied as 
free parameters. In addition, the inner radius of this layer is 
either set to be equal to or varied as a free parameter to 
model a detached C2H2 layer The central star itself is approxi- 
mated with a blackbody of the effective temperature T^ff, which 
is fixed to 2800 K as in the dust shell models discussed above. 
Using the C2H2 opacity given above, the monochromatic inten- 
sity distribution of the C2H2 layer and the star is calculated as 
described in Ohnaka et al. ( I2004bl l. The intensity profile of the 
dust shell (excluding the central star) at each wavelength in the 
band is calculated with the mcsim_mpi code and ray-tracing 
method as in the above dust shell models, and this intensity 
profile is added to that resulting from the central star and the 
C2H2 layer The inner boundary radius of the dust shell is set to 
2.5 Ri,, because the temperature at this radius (1500-1600 K) 
approximately corresponds to the condensation temperature of 
amorphous carbon grains. The spectrally convolved A^-band 
visibiUties are then computed as outlined in Sect.|3] However, it 



has turned out that these models cannot reproduce the shape of 
the 8 jjm absorption feature and the observed A^-band visibili- 
ties simultaneously for the following reason. In these one-layer 
models, the C2H2 layer cannot be too extended, because oth- 
erwise the V4H-V5 and V5 bands would appear in emission, in 
disagreement with the observations. However, with a modest 
geometrical extension of the C2H2 layer, the angular sizes ob- 
served between 8 and 9 //m, which are remarkably larger than 
the photospheric size, cannot be reproduced. 

Therefore, we examine models with a dust shell and two 
C2H2 layers, which are similar to the warm water vapor en- 
velope model for oxygen-rich Miras considered by Ohnaka 
(I2004al l. In these models, the star is surrounded by two lay- 
ers containing C2H2 with different temperatures T^ot and T^ooU 
and the hot and cool layers extend to /?hot and Rcooh respectively 
(see Fig.[TO]for illustration of the model). The temperature and 
the density are assumed to be constant in each of these lay- 
ers. The C2H2 column densities in the hot and cool layers are 
denoted as A^hot and A^cooi, respectively. The computation of in- 
tensity profiles from such two layers are described in Ohnaka 
(I2004bl l. and the total intensity profile is the sum of this com- 
ponent and the dust shell component with its inner boundary 
fixed to 2.5 From these intensity profiles, spectrally con- 
volved visibilities are calculated as in the case of the above 
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Fig. 7. Same as Fig.|6] but for phase 0.49 (minimum light). In the panel a, the filled diamonds and triangles represent the A^-band 
visibilities measured with projected baseline lengths of 124.2 and 41.9 m, respectively, while the solid and dashed lines represent 
the corresponding model visibilities. In the panel b, the solid line represents the model prediction, which is derived by fitting 
the model visibilities predicted for basehnes of 124.2 and 41.9 m with uniform disks. In the panel g, the solid and dashed lines 
almost entirely overlap. 



one-layer model. We search for a best-fit model within the fol- 
lowing parameter range: Thot (K) = 1300 ... 1800 with Arhot 
(K) = 100, Tcooi (K) = 800 ... 1200 K with Ar^ooi (K) = 100, 
Rhot (R*) = 1.1 ... 1.6 with A/Jhot (R*) = 0.1, TJcooi (R*) = 1.2 
... 2.0 with A/?cooi (R*) - 0.1 (with conditions 7?hot < ^cooi). 
logA^hot (cm-2) = 19 ... 22, and logA^cooi (cm^^) = 18 ... 21 
with AlogA^hot = AlogA^cooi = 0.5. The optical depth of the 
dust shell is varied between Ty™ - 0.1 and 0.8 with At™" = 
0.1, while the fraction of SiC, /sic = t^"=/(t™'= + rf^), is var- 
ied from 0.1 to 0.4. 

It should be noted here that not only the wavelength depen- 
dence of the observed visibilities but also the visibility shape 
as a function of spatial frequency constrained by the measure- 
ments at two different baseline lengths are crucial for this mod- 
eling. That is, the effect of the extended dust shell, which is 
over-resolved at the baseline lengths used in our MIDI obser- 
vations, is to lower the total visibility by an amount equal to 
its fractional flux contribution at a given wavelength, as ex- 
plained in Sect. [3] Therefore, the parameters of the dust shell 
roughly defines absolute visibility levels. On the other hand, the 
presence of the C2H2 layers makes the visibility (as a function 
of spatial frequency) steeper at higher spatial frequencies than 
models consisting of a bare central star and a dust shell. The 
C2H2 layers also give rise to absorption or emission at wave- 
lengths between 8 and 9 fim and longward of ~ 1 1 /urn, although 
in this latter wavelength range, the C2H2 feature is blended with 
the SiC emission feature. In this manner, it is possible to con- 
strain the parameters of the hot and cool C2H2 layers as well as 
the dust shell. 



Figures |6l |7] and [8] show the best-fit models consisting of 
C2H2 layers and a dust shell for three epochs observed with 
MIDI, and the parameters of these best-fit models are given in 
Table [3] The figures demonstrate that the models can repro- 
duce the observed A'-band spectra, particularly the absorption 
feature between 8 and 9 jim. At the same time, the observed 
A^-band visibilities — not only the absolute visibility levels but 
also their wavelength dependence — are also fairly reproduced 
by these models. The visibility increase from 8 to lOyum, which 
is seen in all the data sets except for #4 taken with the shortest 
baseline at minimum light, can be well explained by the mod- 
els. The visibilities predicted at 12-13 //m, where the tail of 
the C2H2 V5 band is located, are also in reasonable agreement 
with the MIDI observations. These results suggest that the op- 
tically thick emission (riine ~1-100, see Figs.|6^,|7fe, and[8^) 
from the C2H2 layers, together with the dust emission, is re- 
sponsible for the A^-band angular size significantly larger than 
the photospheric size and that C2H2 plays a role similar to H2O 
in oxygen-rich environments. Also, as will be discussed below, 
this optically thick emission plays an important role in the in- 
terpretation of the controversial apparent depth of the 7 and 
14 jum C2H2 features outlined in Sect.[T] We note that as in the 
case of the dust shell models discussed above, the model vis- 
ibilities are scaled using the angular radius of the central star 
(stellar continuum angular radius) derived by fitting the model 
A^-band spectra to those observed. The stellar photospheric an- 
gular sizes derived in this manner are also given in Table [3] In 
this method, the uncertainty of the stellar diameter is primar- 
ily determined by the uncertainties of the absolutely calibrated 
A^-band spectra of ~10-20%. Based on this, we estimate the 
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Table 3. Derived parameters of the C2H2 (+HCN) layers and 
the dust shell around V Oph. The stellar angular diameter 
(2Ri,/d), where d is the distance of V Oph, was derived by fit- 
ting the model A^-band spectra to those observed. 





phase 0.18 


phase 0.49 


phase 0.65 


Tho, (K) 


1600 + 100 


1600 ± 100 


1400 ± 100 


^hot (^*) 


1.2 ±0.1 


1.2 ±0.1 


1.5 ±0.1 


logTVhotCCjHz) (cm-2) 


20.0 ± 0.5 


20 ± 0.5 


21 ±0.5 


7*0001 (K) 


1000 ± 100 


1000 ± 100 


1000 ± 100 


^cool (R*) 


1.8 ±0.1 


1.4 ±0.1 


1.7 ±0.1 


logA'.oo, (C2H2) (cm-2) 


20.0 ± 0.5 


19 ±0.5 


19.5 ±0.5 


' V 


0.7 ±0.1 


0.2 ± 0.05 


0.5 ± 0.05 


/sic = rS;=/(Tr+T'v"') 


0.2 ±0.1 


0.2 ±0.1 


0.2 ±0.1 


logWho.CHCN) (cm-2) 


< 20 


< 19 


< 21 


logA',„„,(HCN) (cm-2) 


< 20 


< 18 


< 21 


Stellar angular 


4.2 ± 0.3 


4.8 ±0.3 


4.3 ± 0.3 


diameter (mas) 









error of the stellar diameter to be +0.3 mas. It is noteworthy 
that while the observed uniform-disk diameters are the small- 
est at minimum light (see Fig.|2}l), the stellar continuum diam- 
eter does not show such a phase dependence (almost constant 
within the uncertainty or may marginally be the largest at mini- 
mum light). This means that the A^-band apparent size is largely 
governed by the opacity and the geometrical extension of the 
C2H2 layers and the dust shell, and does not reflect the size 
variation of the underlying photosphere. 

The effects of the uncertainty of the effective temperature of 
the central star on the derived parameters are checked by calcu- 
lating models with Teff - 3000 and 2600 K and included in the 
uncertainty of the derived parameters. However, these effects 
have turned out to be relatively minor, because the C2H2 lay- 
ers, particularly the hot layer, are mostly optically thick, mask- 
ing the underlying central star. We also examine the effects of 
HCN on model visibilities and spectra by adding the HCN con- 
tribution to the best-fit models, using the line list of Harris et al 
(I2002i l and the Gaussian line profile with a FWHM of 5 km s ' . 
However, since the spectral feature of HCN is blended with 
the C2H2 V5 band and the SiC dust feature, we cannot put a 
strong constraint on the HCN column densities, and only the 
upper limits of the HCN column densities are derived as given 
in Table [3] 

Figure |9] shows the derived properties of the C2H2 layers 
and the dust shell as functions of phase, and a simple cartoon 
of the models at the three phases is shown in Fig.[TO]to illustrate 
the size of the C2H2 layers and the dust shell and their phase 
dependence. The radius of the cool C2H2 layer show notice- 
able variations: the layer is remarkably more extended at post- 
maximum (phase 0. 18) and post-minimum (phase 0.65) than at 
minimum light (phase 0.49). Furthermore, the column density 
of the cool C2H2 layer appears to be the smallest at minimum 
light, although the uncertainties of the derived column densi- 
ties are rather large. The temperature of the cool layer is found 
to exhibit no temporal variation within the uncertainties of our 
modeling. However, it is not straightforward to interpret these 
temporal variations of the radii as the oscillating motion of the 
C2H2 layer. That is, the smaller radius of the cool layer at min- 



imum light (phase 0.49) does not necessarily mean that C2H2 
gas shrinks from phase 0.18 to 0.49 and bounces back from 
phase 0.49 to 0.65. As discussed in Sect. |5] it is also possible 
that C2H2 layers may be associated with the successive passage 
of shock waves propagating outward and that the observed tem- 
poral variations of the radii and column densities of the C2H2 
layers are simply a series of "snapshots" of such a dynamical 
atmosphere. 

The variation of the C2H2 column densities may appear to 
contradict the phase dependence of the strength of the 8-9 //m 
feature mentioned in Sect. |2| the absorption becomes stronger 
at minimum light (phase 0.49) than at post-maximum (phase 
0.18). This seemingly inconsistent result can be explained as 
follows. At maximum light, the larger geometrical extension 
of the C2H2 layers results in the higher contribution of the 
emission from the extended part of the layers (i.e., beyond the 
central star's limb). In other words, the absorption is partially, 
though not entirely, filled in by this emission component. On 
the other hand, at minimum light, the smaller geometrical ex- 
tension of the C2H2 layers makes the contribution of this ex- 
tended emission component smaller, which results in stronger 
absorption. 

The hot C2H2 layer shows a temporal variation somewhat 
different from that of the cool layer The hot layer becomes 
more extended, denser, and a little cooler at phase 0.65 com- 
pared to the preceding two phases. The dust shell also shows 
noticeable temporal variations: its optical depth is the small- 
est at minimum light, while the fraction of SiC dust does not 
show significant changes within the uncertainty of our model- 
ing. This phase dependence of the optical depth of the dust shell 
is similar to that found in optically bright oxygen-rich Miras 
o Cet and Z Cyg by Suh (I2004I I. who concludes that stronger 
stellar winds near maximum light (higher stellar luminosities) 
lead to higher dust production. We will discuss this point fur- 
ther in Sect.|5] 

We note that while the inner boundary radius of the dust 
shell in the above modeling is fixed by the condensation tem- 
perature of amorphous carbon, the stellar pulsation may lead 
to departure from radiative equilibrium, and dust may form 
closer to the star than it would in radiative equilibrium (Willson 
2000). In order to assess this possibility, we compute dust 
shell models with an inner boundary radius equal to 7?cooi (i-S-, 
no gap between the cool C2H2 layer and the dust shell inner 
boundary) by artificially imposing an upper limit of 1600 K 
on the dust temperatures. We find out that these ad hoc "non- 
equilibrium" models can reproduce the observed A^-band vis- 
ibilities and spectra at three phases if the optical depths are a 
little increased with the parameters of the C2H2 layers kept the 
same as the above models. For example, the C2H2 + dust shell 
model with t™"^ = 1.0 and t^''~ = 0.25 can provide a fit as good 
as the model shown in Fig.|6] This means that the inner bound- 
ary radius of the dust shell cannot be well constrained, because 
the visibility component of the dust shell, which appears as a 
steep drop at the lowest spatial frequencies, is not sampled by 
our MIDI observations using rather long baseline lengths. In 
fact, it would also be possible that dust coexists with the C2H2 
gas. In this case, the innermost part of the dust shell would be at 
least partially masked by the optically thick C2H2 gas, which 
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Fig. 8. Same as Fig.|6l but for phase 0.65 (post-minimum). In the panel a, the filled diamonds and triangles represent the A^-band 
visibilities measured with projected baseline lengths of 87.6 and 120.0 m, respectively, while the solid and dashed lines represent 
the corresponding model visibilities. In the panel b, the solid line represents the model prediction, which is derived by fitting 
the model visibilities predicted for basehnes of 87.6 and 120.0 m with uniform disks. In the panel g, the solid and dashed lines 
almost entirely overlap. 



22 



^.1.5 - 
Si 

m 
3 

K 1 - 



0.5 



0.2 



b. 



1500 - 

£ 

1000 - 



500 



0.2 



0.4 0.6 
phase 



• Hot layer 
oCool layer 

0.8 





21 


'a 










20 


O 








o 


19 



0.2 



a, 

. in 



1.5 - 



• Hot layer 
oCool layer- 



0.4 0.6 
phase 



0.2 



0.4 0.6 
phase 



• Hot layer - 
oCool layer 

0.8 



0.4 0.6 
phase 



5.5 



' 4.5 



' 3.5 



0.2 



0.4 0.6 
phase 



0.8 



Fig. 9. Temporal variations of the radii, temperatures, and column densities of the hot and cool C2H2 layers (a, b, and c), the 
optical depth of amorphous carbon dust in the dust shell (d), and the continuum stellar angular diameter (e). 



would mimic the dust shell with its inner boundary equal to 
Rcooi- It is also worth noting that the non-LTE grain formation 
calculation of Kozasa et al. (1 19961 1 shows that SiC grains form 
first and then amorphous carbon mantles form around the SiC 
cores at low mass loss rates, which is not considered in our 



present model. These dust shell properties could be better con- 
strained by observations with much shorter baselines (^20 m). 

There are still discrepancies between the best-fit models 
and the observed data. At phase 0.18 and 0.65, the predicted 
visibilities (as a function spatial frequency) are still somewhat 
too flat compared to the MIDI visibilities measured at two base- 
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Fig. 10. A simple cartoon illustrating the best-fit models at phases 0.18, 0.49, and 0.65. The sizes of the components at three 
phases are to scale except for the outer radius of the dust shell. 



line lengths. This means that the C2H2 layers may be larger 
than in these models. In the framework of our two-layer mod- 
eling, however, adoption of larger radii for the hot and/or cool 
C2H2 layers turns the absorption feature between 8 and 9 /im 
to emission. Moreover, while we assume for simplicity that the 
hot C2H2 layer extends down to the photosphere and the hot 
and cool C2H2 layers are in contact, it is possible that these 
C2H2 layers are detached and geometrically thin. This may also 
explain the discrepancy between the models and the observed 
data. In any case, our two-layer model with constant tempera- 
tures and densities is a rather simplistic representation of real 
density and temperature distributions, and more sophisticated 
models of the C2H2 layers constrained by further complimen- 
tary observational data may reproduce the MIDI data better 
For example, mid-infrared high-resolution spectra of the C2H2 
V4-1-V5 and V5 bands would shed more light on the physical 
properties of the C2H2 layers, as 11 jum high-resolution spectra 
have proven to be strong observational constraints in the case 
of oxygen-rich Mira variables (e.g., Ohnaka 2004a; Ohnaka 
12006b J . Also, although we assume LTE in the calculation of the 
molecular opacities, it is necessary to examine non-LTE effects 
when an extensive high-temperature line list of C2H2 becomes 
available. 

5. Discussion 

Yamamura et al. (119971) analyzed the mid-infrared spectra of 
several carbon stars obtained with ISO/SWS and derived the 
C2H2 column densities and temperatures from the 13.7 jjm fea- 
ture. The C2H2 column densities of lo"^"^*^^ cm"^ are somewhat 
lower than those we derived for the cool C2H2 layer, while the 
temperatures of 1200+300 K that they derived are in fair agree- 
ment with those of the cool layer of our models. It should be 
noted that they measured the strength of the 13.7 jum feature 
with respect to the local continuum levels, which are affected 
by the broad (hot) C2H2 absorption component as well as dust 
emission. Therefore, the C2H2 column density and temperature 
derived by Yamamura et al. (1 19971 1 may refer to the cooler, out- 
ermost layers. 

Recently, Nowotny et al. (I2005at l2005bl) have presented a 
self-consistent dynamical model including dust formation pro- 
cesses for the carbon-rich Mira S Cep (M k 10"^ M© yr"'). 
They demonstrate that their model can qualitatively reproduce 
the observed temporal variations of the radial velocity and line 
profiles of infrared molecular lines such as CO and CN. The 



density distribution of their model shows a number of step-like 
features resulting from the formation of shock fronts propagat- 
ing outward. It is plausible that polyatomic molecules such as 
C2H2 and HCN form behind shock fronts — possibly at the in- 
nermost ones where density is high — just as H2O shells form 
in oxygen-rich Miras (Tej et al. l2003l ). Therefore, it is well pos- 
sible that the observed temporal variation of the radii of the 
C2H2 layers results from changes of the positions of C2H2 for- 
mation behind shock fronts. This means that the cool layer ex- 
tending to ~1.8 Ri, at phase 0.18 moves farther outward and 
may no longer contribute to the A^-band flux at phase 0.49. And 
the small radius of the cool layer obtained at phase 0.49 sug- 
gests that a new C2H2 formation front may just be emerging 
at a smaller radius and that this front may build up and move 
outward from phase 0.49 to 0.65. 

The phase dependence of the grain optical depths we de- 
rived for V Oph can also be interpreted in a similar manner. The 
model presented by Nowotny et al. (I2005at I2005bl l show that 
dust is produced in discrete shells behind shock fronts where 
the density is sufficiently high and the temperature is below 
condensation temperature. If such a hot dust layer is respon- 
sible for the higher grain optical depths we derived for phase 
0.18 (post-maximum light), the dust layer is expected to expand 
from phase 0.18 to 0.49 (minimum light) and become diluted. 
At minimum light, a new (inner) dust layer can be forming, but 
still with a low degree of condensation, which can qualitatively 
explain the smallest grain optical depths at minimum light. 
When this new dust layer has fully grown, it can provide higher 
grain optical depths at post-minimum phases. However, as ex- 
plained above, our present MIDI data cannot constrain the in- 
ner boundary radius of the dust shell strongly enough to detect 
the apparent motion of dust layers expected from this scenario. 
We also note that the mass loss rate of V Oph is of the order of 
10"^ Mq yr"', which is by two orders of magnitude lower than 
that of S Cep. It is necessary to compare with dynamical mod- 
els computed with parameters corresponding to V Oph for bet- 
ter understanding physical processes responsible for molecule 
and dust formation close to the star: whether the extended 
C2H2 layers can be explained by dust-driven winds triggered 
by large-amplitude stellar pulsation alone or some other mech- 
anism is operating such as Alfven wave-driven winds (e.g., 
Vidotto & Janteco-Pereira l2006l Suzuki |2006] |. 

The result that our two-layer models can reproduce the ob- 
served C2H2 absorption feature between 8 and 9 jim implies 
that the emission from these dense, geometrically extended 
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C2H2 layers may be responsible for the observed weakness of 
the C2H2 and HCN absorption features at 7 and 14 /im. While 
the dust thermal emission also contributes to rendering these 
molecular features weak, the low mass loss rate of V Oph sug- 
gests that the optically thick emission from the C2H2 layers 
plays a major role at least in optically bright Mira-type car- 
bon stars. In order to confirm this, it would be crucial to com- 
pare model predictions with observed strengths of the C2H2 
and HCN features and their temporal variations at 3 fim (ac- 
cessible from the ground) as well as 7 and 14 foa (with Spitzer 
and/or AKARI). 

6. Concluding remarks 

We have carried out multi-epoch A^-band spectro- 
interferometric observations of the carbon-rich Mira V Oph 
and spatially resolved the object between 8 and 13 yum. The 
observed A^-band uniform-disk diameters are by a factor of 
2-3 larger than the estimated stellar continuum size, as has 
already been established for oxygen-rich Miras. The visibilities 
obtained at phases 0.18, 0.49, and 0.65 reveal that the object 
appears smaller at phase 0.49 (minimum light) than at phases 
0.18 and 0.65. Our simple modeling shows that the observed 
A^-band visibilities and spectra can be reproduced by optically 
thick C2H2 layers and a dust shell consisting of amorphous 
carbon and SiC. The cool C2H2 layer is more extended and 
denser at phases 0.18 and 0.65 than at 0.49, while the hot 
C2H2 layer becomes more extended and denser at phase 0.65 
with no noticeable change between phase 0.18 and 0.49. This 
modeling also implies that the weakness of the C2H2 and HCN 
absorption features at 7 and 14 yum observed in carbon-rich 
Miras can be explained by emission from these extended, 
dense C2H2 (and HCN) layers, together with dust emission. 
Since our modeling is based on the C2H2 opacity computed 
with the RRHO band-model approximation, it is essential to 
examine the present results with an extensive high-temperature 
line list of C2H2 in future. However, while the C2H2 column 
densities can be directly affected by the uncertainty of the 
C2H2 opacity, the trend of the phase-dependent variation of the 
C2H2 column densities found in the present work is less likely 
to be largely affected. 

In order to better understand the phase and cycle depen- 
dence of the physical properties of the outer atmosphere and the 
dust shell, it is indispensable to carry out near-infrared interfer- 
ometric observations simultaneously with MIDI observations. 
Such coordinated observations are in fact feasible with VLTI, 
using MIDI and the near-infrared beam combiner, AMBER, 
which is capable of measuring visibilities in the J, H, and K 
bands with three different spectral resolutions of 75, 1500, and 
12000. With the highest spectral resolution, it is possible to 
measure the stellar continuum diameter, although not entirely 
free from molecular line contamination, and to obtain a more 
comprehensive picture of how the outer atmosphere and the 
dust shell respond to stellar pulsation. 
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